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Abstract 

The purpose of this paper is to expound and clarify the mathematics 
and explanations commonly employed in certain notable areas of astron- 
omy and astrophysics. The first section concentrates upon the mathemat- 
ics employed to represent and understand stellar structure and evolution. 
The second section analyses two different explanations for the structure 
of spiral galaxies. 

1 Stellar structure and evolution 

A star forms when a body of gas contracts under its own gravitational attraction, 
and the pressure and temperature created at the centre of the agglomerated mass 
is sufficient to ignite the nuclear fusion of atomic nuclei. The radiation released 
as a by-product of the nuclear fusion eventually reaches the surface of the star 
and becomes starlight. 

Newtonian astrophysics represents the gaseous content of a star as a fluid 
which occupies a compact open subset $7 C M'^ of three-dimensional Euclidean 
space. The fluid is considered to possess a mass density scalar field p and a 
pressure scalar field p. To study stellar structure in Newtonian astrophysics, it 
is common to assume that the star is spherically symmetric and static (time- 
independent). The requirement that the star be static is equivalent to the 
requirement that the star is in hydrostatic and thermal equilibrium. For a 
static spherically symmetric star, the fluid occupies a soHd ball of radius r, the 
mass density and pressure are time independent, and the velocity vector of the 
fluid vanishes. Using spherical polar coordinates, (r, 0, 4>), the pressure and the 
density are independent of the angular coordinates (0, (/>), but vary as a function 
of the radius r. The fluid is assumed to be self-gravitating and to satisfy the 
Poisson equation of Newtonian gravitation, 

= -47rGp , 

with respect to a gravitational potential scalar field $. is the Laplacian, also 
denoted as A in many texts. 
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There arc four differential equations which govern the structure of a static, 
spherically symmetric star. The four equations taken together constitute a cou- 
pled ('simultaneous') set of ordinary first-order, non-linear differential equa- 
tions. We now proceed to introduce and discuss these equations. 

Assuming spherical symmetry and time-independence, the mass density is 
a function of radius alone p{r), and from the definition of mass density and 
Euclidean geometry, it follows that m(r), the mass enclosed within the surface 
of radius r, is given by the expression: 



m(r) = / 4Tr{r'fp{r')dr' 
Jo 



(The prime here simply indicates the use of a dummy variable for the integra- 
tion.) This gives us our first differential equation for stellar structure, sometimes 
called the mass continuity equation: 

dm(r) ^ 2 / \ 
dr 

The Newtonian gravitational potential for a self-gravitating, spherically sym- 
metric body satisfying the Poisson equation, is 4> = —Gm{r)r~^ . Given a grav- 
itational potential the gravitational force exerted per unit mass is specified 
by the vector field — grad and the gravitational force per unit volume is spec- 
ified by F = — p grad Most texts denote the gradient operator as V, but this 
can be confused with the covariant derivative operator in differential geometry. 
Whilst the gradient of a scalar field is a contravariant vector field, the covariant 
derivative of a scalar field is a covariant vector field. Hence, to avoid confusion, 
we shall write 'grad' rather than V. Given a choice of coordinates {x\, ...,a;„) 
on a manifold, and given a metric tensor field gijdx^ (g) dx^ which specifies the 
geometry of space, the gradient of an arbitrary scalar field / can be expressed 
as 

grad/ = 5'^|^T^, 

oxi oxj 

where g^^ is the inverse of the matrix gij . In the spherical polar coordinates we 
have chosen, {xi,X2, X3) = (r, 6, </>), the flat Euclidean metric takes the form 

d,s^ = dr (g) dr + r^d0 (E) dO + sin^ 9d(j) (g) d(j) , 
hence the gradient of / takes the form 

Assuming spherical symmetry of the gravitational potential $ entails that 
d^/d9 and d'^/dcj) vanish, hence 



2 



1 * d'^ d 
grad $ = — — 
dr or 

_ d{-Gni{r)r~^) d 

dr dr 

— Gm(r)r^^ — . 

or 

The requirement of hydrostatic equihbrium corresponds to the requirement 
that the inward gravitational force —p grad $ is balanced at each point by the 
outward pressure gradient force —grad 

grad p + p grad $ = . 
Spherical symmetry entails that dp/ 86 and dp/dcf) vanish, hence 

, dp d , \ , \ -2 

gradp = -r-^ — ~Gm(r)p(rjr — . 

dr or Or 

This gives us our second differential equation for stellar structure: 

dP{r) Gm{r)p{r) 
dr 

Let l(r) denote the amount of energy passing through the surface of radius 
r per unit time. Let e(r) denote the energy production coefficient, the amount 
of energy released at radius r in the star per unit time and per unit mass. It 
follows that 47rr^p(r)e(r) is the energy released per unit time at radius r. A star 
in thermal equilibrium is an open system for which there are large energy flows 
out of the system, and in which there are significant temperature and pressure 
gradients, but in which, nevertheless, the temperature and pressure profile of 
the star remain constant. For a star in thermal equilibrium, none of the energy 
released is used to heat up the star or change its volume. Hence, for such a star, 
the rate of energy flow at radius r is given by the expression: 

l{r) ^ [ 4:7r{r''fp{r')e(r')dr' . 
Jo 

This gives us our third differential equation of stellar structure, 

^^iM ^ Anr^ p{r)e{r) . 
dr 

For a star of radius R the luminosity is defined to be L = 1{R), the amount of 
energy passing through the outer surface of the star per unit time. When a star 
resides in a state of thermal equilibrium, the total amount of energy produced 
by the star per unit time equals the amount of energy radiated from the outer 
surface per unit time. Hence, the luminosity L of such a star is given by the 
equation: 



^The pressure at each radius is a sum of the gas pressure and the radiation pressure. 
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R 

2 , 



L = 47rr p(r)e{r)dr . 
Jo 

The energy per unit time l{r) passing through the sphere of radius r can be 
expressed as 

l{r) = 4TTr^F{r) , 

where F is the energy flux, the energy flow per unit time per unit area. When 
there is a temperature gradient in a body, energy flow will occur by conduction 
and radiation at the very least. Assuming, for simplicity, the absence of con- 
vection, the energy flux in a star can be broken down into a conductive energy 
flux Fcondif) and a radiative flux Frad{f)- In both cases, the energy flow is 
proportional to the temperature gradient as follows (Tayler 1994, p63): 



F,,^Jr) = 4acr3(r) dT{r) 
5Kcond{r)p{r) dr 



and 



Fradir) 



iacT^ir) dT{r) 



^Krad{r)p{r) dr 

Kcond, the coefficient of conductive opacity, measures the resistance to the flow of 
heat by conduction, and Krad, the coefficient of radiative opacity, measures the 
resistance to the flow of heat by radiation, a is the radiation density constant. 



a = 7.57 X 10~^^ erg cm~^K 



Assuming that energy flow in a star is due to conduction and radiation^ one 
obtains (Tayler 1994, p64): 

9. .X , XX Wiracr'^T^ir) dT(r) 



3K(r)p{r) dr 



where 

1 1 1 



li{r) Kcondir) Kradir) 

Re-arranging, one obtains the following differential equation for the temperature 
gradient in a star, which is the fourth equation governing stellar structure: 

dT{r) 3K(r)p(r) l{r) 



dr AacT^{r) Anr'^ 

The chemical composition of a spherically symmetric and static star can be 
specified by a set of functions X.i{r), which specify the fraction of unit mass 
consisting of nuclei of type z, for i = 1, ...,/. As such, EiXi(r) = lH 

^Conduction is only significant in 'compact' stars, i.e., white dwarf and neutron stars. 

•^It is a common approximation to use only three fractions {X, Y, Z), which are such that 
X + Y + Z = I, and which specify, respectively, the proportion of hydrogen, the proportion 
of helium, and the proportion of all other elements in the constitution of a star. 
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Before the equations of stellar structure can be solved, one must specify an 
equation of state for the pressure, 




e = e(p,T,Xi, 




and an expression for the opacity 



These expressions are sometimes called the constitutive equations. 

If these three functions are specified, and four boundary conditions are spec- 
ified, then one can solve the four differential equations of stellar structure. The 
four boundary conditions can be chosen to be: 



Note here that the star is idealised as an open solid ball of radius R, hence the 
boundary of the star at radius R is not itself considered to be part of the star. 
Choosing p = (or P = 0) and T = at a finite radius r = i? is merely a simple 
idealisation. In reality, the density, pressure and temperature of a star gradually 
descend to the non-zero value of the interstellar medium in the neighbourhood 
of the star. 

Solving the four differential equations for stellar structure then gives the 

mass, pressure, energy flow and temperature as a function of radius. However, 
in practical terms, one might wish to specify the total mass M and chemical 
composition of a star, and then obtain, amongst other things, the radius of 
such a star in hydrostatic and thermal equilibrium. It is therefore conventional 
to re-cast the differential equations with mass m, rather than radius r, as the 
independent variable. Doing so, one obtains (Tayler 1994, p70-71): 



m = 0, / = atr = 



and 



p = 0, T = at r = R. 



dr{m) 



1 



dP{m) 



dm 



A-Kr'^{m)p{m) ' 
Gm 



dm 



47rr^(m) ' 



dl{m) 



= e(TO) , 



dm 



dT{m) 



ZK.{m)l{m) 



dm 



647r2acr4(m)r3(m) ' 
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The chemical composition of the star is then expressed in terms of functions 
{X,{m):i = l,...,I}. 

One can then choose boundary conditions for the re-formulated equations of 
stellar structure, such as: 

r = 0, 1 = atTO = 0, 

and 

p = , T = Q at TO = M . 

Again, one could choose P = at to = M rather than p = 0, and, again, this is 
merely the simplest idealisation on offer. 

Now consider the relationship between theory and observation. The ob- 
servational state of a star at a moment in time can be specified by just two 
parameters: its luminosity and 'surface' temperature. Let us take these in turn: 

1. If one calculates the distance to an observable star, say by parallax, then 
one can infer the luminosity of the star from its apparent luminosity. 

2. Whilst the temperature of the surface of a star was idealised in the bound- 
ary conditions above to be zero, a more sophisticated model divides a star 
into its interior and atmosphere. The lowest level of the atmosphere is 
called the photosphere, and the 'surface' temperature of a star is deemed 
to be the temperature of the photosphere. The temperature of a star's 
surface, in this sense, is intimately related to the spectral type of the star, 
defined by the absorption lines in the spectrum of light emitted by the 
photosphere. Given that the temperature of a star's photosphere largely 
determines the type of its spectrum, the observational state of a star can 
be specified by its luminosity and spectral type. There are ten spectral 
types which, in order of decreasing temperature, are referred to as O, B, 
A, F, G, K, M, R, N, and S. 

If one knows both the luminosity and temperature of a star, then one can cal- 
culate the radius of the star. The luminosity of a star is a function of the radius 
R of the star and the effectiv^ surface temperature Tg of the star, according to 
the following expression, 

L = AirR^aT^ , 

where a is the Stefan-Boltzmann constant, related to the radiation density con- 
stant a by CT = ac/4. Given L and Tg, one can obviously use this expression to 
calculate R. 

The observational state of a star can be represented by a point on the 
Hertzsprung-Russell diagram, a two-dimensional rectangle, coordinatized by lu- 
minosity on the vertical axis, and either spectral type or temperature on the 

* Strictly, the effective surface temperature Te is defined as the temperature of a black body 
with the same radius and luminosity as the star. 
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horizontal axis. The observational history of a star traces out a path on the 
Hertzsprung-Russell diagram. 

The luminosity and temperature of a star in hydrostatic and thermal equi- 
librium are the observational properties which can be explained by solving the 
differential equations of stellar structure. This explanation falls under the aegis 
of the deductive-nomological (D-N) account of scientific explanation. In such 
explanations, one explains certain phenomenal facts by logically deriving them 
from the conjunction of general laws and particular specified circumstances. In 
this case, the mass and chemical composition of a star are the specified circum- 
stances, and given a specification of the constitutive equations, one can explain 
the luminosity and temperature of a star from the conjunction of the equations 
of stellar structure and the given mass and chemical composition. 

The initial mass and chemical composition of different stars vary, and these 
are the characteristics which determine the history and lifetime of a star on the 
Hertzsprung-Russell diagram. Stars which are burning hydrogen into helium, 
in a state of hydrostatic and thermal equilibrium, occupy a roughly diagonal 
channel, running from the top left to the bottom right of the Hertzsprung- 
Russell diagram, called the main sequence. The initial mass, and to a lesser 
extent, the initial chemical composition of a star, determine its initial location 
in the main sequence. Stars with greater initial mass are more luminous, and 
occupy positions further up the main sequence slope. The larger the mass and 
luminosity of a star, the shorter the length of time it will spend on the main 
sequence until its hydrogen fuel is expended. 

The Vogt-Russell 'theorem', conjectured independently by Heinrich Vogt 
(1926) and Henry Norris Russell (Russell et al 1927, p910), states that with the 
total mass and chemical composition specified, the equations of stellar struc- 
ture admit a unique solution. If true, the Vogt-Russell theorem would entail 
that each mass and chemical composition corresponds to a unique equilibrium 
configuration. This, of course, is consistent with the fact that stars of the same 
mass but different chemical profile, can possess different equilibrium states. The 
chemical composition of a star will change during its lifetime as nuclear fusion 
converts lighter elements into heavier elements, and a star will tend to pass 
through a sequence of different equilibrium configurations. 

However, whilst the Vogt-Russell theorem is essentially true, it is not strictly 
true, and it is certainly no theorem. A system of four coupled, non-linear, ordi- 
nary differential equations, with boundary conditions specified at two different 
points, does not necessarily admit a unique solution. Whilst there is usually a 
unique equilibrium configuration for a given mass and chemical composition, for 
some combinations of mass and chemical composition there are multiple solu- 
tions to the equations of stellar structure. The notion that there was actually a 
'theorem' was promulgated in textbooks, and never given a rigorous proof. As a 
conjecture, it is now known to be false. In fact, both the existence and unique- 
ness parts of the conjecture fail. For example, there is no solution composed 
of helium with a mass less than about 0.3 solar masses, and Cox and Salpeter 
(1964) demonstrated that there are two different solutions for stars of the same 
mass, burning purely helium, close to this minimum mass. This is referred to 
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as the double-valued helium main sequence, in the sense that, for some range of 
masses, two values of radius R are possible for each value of mass. Kahler (1978) 
tracked the post main sequence evolution of the helium core of a star of two 
solar masses, and found different equilibrium configurations of different radius, 
but with the same mass. Gabriel and Noels-Grotsch (1968) demonstrated that 
the minimum mass of a solution composed entirely of carbon is about 0.9 solar 
masses, and that, once again, there are two possible solutions for each mass 
value greater than the minimum mass. In one branch of such double- valued 
main sequences, the radius and luminosity increase with increasing mass, whilst 
in the other branch, there are smaller radii, and the luminosity decreases with in- 
creasing mass (Hansen 1978, p23). The latter, anomalous branches, correspond 
physically to the presence of electron degeneracy. It is therefore not strictly 
true to say that the mass and chemical composition of a star in equilibrium 
uniquely determine its structure. However, as Hansen et al point out, "the idea 
of uniqueness is still useful in that among a set of models all having the same 
mass and run of composition, usually only one seems to correspond to a real 
star or to have come from some realistic line of stellar evolution. The others are 
unstable in some fundamental way (as far as we know)," (2004, p331). 

Given a fixed mass and chemical composition, each solution to the differ- 
ential equations of stellar structure corresponds to a pair of values (P^T^) = 
(P(0),T(0)) for the central pressure and temperature of the star, or, equiva- 
lently, to a pair of values {B.,L) = {r(M), L(M)) for the surface radius and 
luminosity. The significance of this is that the radius R and luminosity L of a 
star are both inferrable from observation. If the mass and chemical composi- 
tion determined a imique solution, then the mass and chemical composition of a 
star in equilibrium would uniquely determine its radius and luminosity {R,L). 
If, on the contrary, there is no unique solution for a particular combination 
of mass and chemical composition, then there are multiple corresponding pairs 
{{R,L)i : i = l,...,n}. 

If one drops the requirements of hydrostatic and thermal equilibrium, then 
one can obtain the equations of stellar structure and evolution which govern the 
history of a Newtonian star. Assuming spherical symmetry, the equations are 
(Kippenhahn and Weigert 1990, p64): 

dr{m, t) 1 



dm A'Kr'^{m,t)p{m,t) ' 

dP{m,t) _ Gm 1 d'^r{m,t) 

dm ~ ~A'Kr^{m,t) ~ 4m^n^t) dt^ ' 

dl(m,t) _ dT{m,t) S{m,t) dP{m,t) 

dm ^ ^ dt p{m,t) dt ' 

dT{m,t) 2>K{m,t)l(m,t) 
dm ~ ~647r2ocr4(m,t)T3(m,f) ' 

Notice that these are partial differential equations, whilst the equations of stellar 
structure are ordinary differential equations, cp is the specific heat at constant 



8 



pressure and S = —{dlnp/dlnT)p, (see Kippcnhahn and Weigert, pl9). The 
fourth equation here continues to assume that energy transport is due to radi- 
ation and conduction alone, but this equation can be generaUzed. 

Timc-dcpcndcncc requires one to introduce additional equations for the time- 
dependence of the mass fractions Xf. 

dX,{m,t) rrij 

TOj is the mass of the nuclei of type i, and is the rate at which nuclei of type 
i are transformed into nuclei of type j per unit volume. 

Solving these equations of stellar structure and evolution requires the spec- 
ification of initial conditions as well as boundary conditions. 

The histories of different star types can intersect when plotted as paths 
on the Hertzsprung-Russell diagram. For example, after a star of around one 
solar mass has expended all the hydrogen in its core, it leaves the main sequence, 
increasing in radius and decreasing in temperature, to the effect of a net increase 
in luminosity. This evolutionary path takes the star up the Hertzsprung-Russell 
diagram to become a red giant. It then subsequently increases in temperature, 
crossing the main sequence itself from right to left, before it sheds its outer 
layers in the form of a so-called planetary nebula. Where this path crosses the 
main sequence, the mass of the star is considerably below the mass of stars at 
that point on the main sequence. If two different types of star can occupy the 
same point on the Hertzsprung-Russell diagram, it entails that knowledge of 
the luminosity and temperature of a star alone is not sufficient to determine the 
unique evolutionary history of a star. The Hertzsprung-Russell diagram docs 
not provide a state space for stars in the sense of a state space ('phase space') 
in Hamiltonian mechanics. 

Nuclear fusion inside stars is responsible for creating almost all of the atoms 
in our universe which are heavier than hydrogen or helium. These heavier 
elements are gencrically called 'metals' by astrophysicists. During their lifetime, 
and particularly at the end of their lifetime, stars will eject a proportion of 
their mass into the interstellar medium, and this mass will contain, in some 
proportion, the metallic elements created by fusion. The mass returned into the 
interstellar medium is re-cycled in the formation of subsequent generations of 
stars. 

The stars in our galaxy are classified, in a coarse-grained fashion, as either 
Population I or Population II. The Population II stars were the first generation 
of stars formed in our galaxy, and as such, are metal-poor stars. Population 
I stars were formed from an interstellar medium which already contained the 
metallic elements created by the first generation of stars, hence the Popula- 
tion I stars have a greater proportion of metallic elements in their chemical 
constitution. 
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Figure 1: Spiral galaxy M74. 



2 Galaxies 

A galaxy is a collection of ^ 10^^ stars, which is gravitationally bound together, 
and which is ^ 10^ light years in diameter. Stars within a galaxy are separated 
by distances of ~ 10 light years, (Rindler 2001, p350). There are two main 
types of galaxy: spiral and elliptical. Most of the stars in a spiral galaxy are 
concentrated in a flattened disk, which rotates about its center. The rotation, 
however, is not rigid, and the angular velocity of the rotation decreases as a 
function of radius. The average rotation period in a spiral galaxy is about 
100 million years, (ibid.). The disk of a spiral galaxy is surrounded by a more 
diffuse 'halo' of stars, distributed in a spherically symmetric fashion. The halo 
is populated only by Population II stars. Whilst successive generations of stars 
form throughout most of the disk in a spiral galaxy, there is typically a bulge at 
the center of the disk, and this bulge also tends to contain Population II stars. 
At the very center of the disk there is typically a supermassive black hole. A 
spiral galaxy possesses an interstellar medium which consists of gas and dustH 
The stars in the disk of a spiral galaxy are not actually distributed in spiral 

^In relativistic cosmology, 'dust' means a pressure-less fluid. In astronomy and astro- 
physics, however, dust means tiny grains of solid matter (Nicolson 1999, pl57). 
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patterns. Rather, the spirals trace the brightest regions in the disk of the galaxy. 
The brightest regions are those which contain both the O- and B-type stars and 
the HII regions of ionized hydrogen, created by the ultraviolet radiation from 
those high-temperature O and B-type stars. Because the brightest (high mass) 
stars are also the stars with the shortest lifetimes, the bright high-mass stars 
can only be found in regions of recent star formation, and these regions happen 
to be the spiral arms. 

In an elliptical galaxy, the stars are uniformly distributed within a 3- 
dimensional ellipsoid, and there is no overall rotation. In fact, the velocities 
of the stars in an elliptical galaxy are randomly distributed in all directions. 
Much like the halo of a spiral galaxy, there appears to be only one generation 
of star formation in an elliptical galaxy. 

Our own galaxy, the Milky Way, is a so-called barred-spiral, in which the 
spiral arms emanate from the ends of an apparently rigid bar of luminous matter, 
(Nicolson 1999, p206). There appear to be two main spiral arms in the Milky 
Way: the Norma arm and the Sagittarius arm. In addition, there are a number 
of smaller arms, and the Sun can be found within the Orion arm or 'spur', which 
lies outside the Sagittarius arm, and inside the Perseus arm (Nicolson p202). 
The Sun resides inside a 'bubble' three-hundred light years in diameter, in which 
the interstellar medium has been cleared by a supernova explosion (Smolin 1997, 
pl24). Prom this perspective, of the 10^^ stars in our galaxy, only about 7000 
can be seen with the naked eye, (Rindler p35G). Moreover, only 1% of the light 
which falls upon the Earth comes from beyond our galaxy, (Disney 2000, p4). 

The mutual gravitation of the stars in a galaxy define an escape velocity, and, 
at any one time, most of the stars in a galaxy will be moving with a velocity 
less than the escape velocity. In this sense alone, a galaxy cannot be treated 
as a collection of isolated, independent entities. Rather, a galaxy must be 
represented as a system, in the sense that it constitutes a collection of mutually 
interacting parts. In a spiral galaxy in particular, one can treat the stars as the 
members of a population, and the interstellar medium as the environment with 
which that population interacts, via feedback loops. In this sense, it has been 
suggested that a spiral galaxy can be treated as an ecological system, (Smolin 
1997, Chapter 9). Perhaps, however, it would be more accurate to say that 
biological ecosystems and spiral galaxies arc both instances of the same type of 
formal system. In other words, there is a type of formal system which contains a 
population, an environment, and a set of feedback loops between the population 
and the environment, and spiral galaxies provide instances of this system- type 
just as much as biological ecosystems do. 

Smolin asserts that "There are processes by which the matter of the inter- 
stellar medium is converted into stars and there are processes by which matter 
is returned from the stars to the interstellar medium. To understand what a 
galaxy is, and especially to understand it as a system, is then primarily to un- 
derstand the processes that govern the flow of matter and energy between the 
stars and the interstellar medium," (1997 pll8). He points out (pl24) that the 
interstellar medium is a system far from thermodynamic equilibrium, consisting 
of a number of different components of different densities, temperatures, and 
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compositions. The different components include normal atomic gas; cold, dense 
Giant Molecular Clouds (GMCs); and regions of hot, dilute plasma, otherwise 
known as HII regions, denoting ionized hydrogen. The relative amount of mat- 
ter in these different components remains approximately constant over time. 
Smolin asserts that for a system to be in such a stable, but far-from-equilibrium 
state, there must be processes which cycle the material among the different 
components, and the rates of these processes must be controlled by feedback 
mechanisms. 

The arms of a spiral galaxy trail behind with respect to the direction of 
galactic rotation, hence the obvious explanation for their shape is simply the 
fact that the angular rotation velocity decreases with radial distance from the 
centre of a galaxy. However, at the very least, this differential rotation cannot 
be the sole explanation because spiral galaxies have existed for at least 10 billion 
years, and given a typical rotation period of 100 million years, this means that a 
typical spiral galaxy has completed at least 100 rotations, and "would long ago 
have wrapped its arms to extinction if they had been created early in its history 
and had not been sustained by some ongoing process," (Nicolson 1999, p208). 
Observations also suggest that the spiral arms do not rotate with the galaxy, but 
at a slightly slower rate. Note carefully, however, that an explanation which uses 
differential rotation in part is not excluded; it may still be possible to explain 
the spiral arms by invoking differential rotation and an ongoing process which 
continually creates the arms. What cannot be done is to explain spiral arms 
simply in terms of differential rotation acting upon the initial conditions under 
which galaxies are created. The burden of explanation has been shifted from the 
initial conditions under which a spiral galaxy is created to the ongoing processes 
which continually operate within such a galaxy. 

There are two processes which have been suggested to explain the ongoing 
formation of spiral arms: density-waves and self-propagating star formation. 
As Nicolson (p208) points out, both processes may participate in spiral arm 
formation, and may be of differing importance in different types of spiral galaxy. 
There are so-called 'grand-design' spirals, which possess thin, long and well- 
defined spiral arms (see Figure [2); there are flocculent spirals, which possess 
many fluffy, poorly-defined spiral arms (see Figure [3]); and there is a continuum 
of intermediate cases, of which the Milky Way is one such. Karttunen et al 
point out that "in multiarmed galaxies the spiral arms may be short-lived, 
constantly forming and disappearing, but extensive, regular, two-armed patterns 
have to be more long-lived," (2003 p359). As Nicolson suggests, density- wave 
mechanisms are more appropriate for grand-design spirals, and self-propagating 
star formation is more appropriate for flocculent spirals. 

It was Bertil Lindblad who suggested in the 1920s that spiral arms could 
be produced by density waves propagating in a galaxy. The peaks and troughs 
of a density wave are independent of any particular elements from the medium 
through which the density wave is propagating. Thus, the density- wave expla- 
nation accepts that stars pass in and out of the spiral arms. However, it also 
accepts that spiral arms are genuine indicators of higher-than-average stellar 
density. The postulated density wave travels faster than the speed of sound in 
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Figure 2: Hubble Space Telescope image of M51, the Whirlpool galaxy. 



the interstellar medium, and the shock wave purportedly triggers compression of 
the interstellar medium, and the observed star formation in the spiral arms. The 
spiral arms therefore trace the brightest regions in the disk of a galaxy because 
the brightest regions are those which contain the shortlived O- and B-type stars. 
As Tayler states, "the most massive stars have such a short lifetime that, by the 
time they have ceased to be luminous, the spiral pattern has hardly changed its 
position. They should therefore only be found in the spiral regions. . . Stars of 
lower mass, such as the Sun, have a main sequence lifetime which is equal to 
many rotation periods of the pattern. Such stars should therefore be observed 
throughout the disk with no significant correlation with the present position of 
the spiral arms," (1993, pl45). 

C.C.Lin and Frank Shu developed the density-wave idea further in 1964 by 
postulating that the stars in a spiral galaxy travel in slightly elliptical orbits 
which precess with the passage of timeO This postulate was complemented by 
the suggestion of J. Kalnajs in 1973 that if the orientation of the major axis of 
these ellipses varies by a small angular increment at increasing distance from 
the galactic centre, then the ellipses fail to be uniformly separated, and where 

^See Pasha (2004) for a detailed history of these ideas. 
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they come close together, they produce the appearance of grand-design spiral 
arms. 

Kalnajs's explanation of spiral arms can be seen as a structural explanation, 
in the sense that the spiral arms are represented as structural elements in a 
geometrical model, and these structural elements are linked to other structural 
elements in the model, namely the precessing elliptical orbits, each aligned at a 
small angle to the one inside it. 




Figure 3: Flocculent spiral galaxy NGC4414. 



Nicolson suggests that density-wave spiral arms "are probably sustained ei- 
ther by the asymmetric gravitational field associated with a central bar struc- 
ture (typical of grand-design spirals) or by gravitational disturbances caused by 
neighboring galaxies, or by a combination of both," (1999, p208). 

In the 1970s Mueller and Arnett suggested an alternative explanation for 
spiral arms, which rejects the idea that all spiral arms are density waves, and 
postulates instead that at least some are self-propagating waves of star for- 
mation. Unlike density waves, these waves purportedly do not require one to 
postulate external mechanisms, such as the asymmetric gravitational field of the 
central bar, or the tidal forces caused by other galaxies. Let us, then, examine 
star formation processes in a little more detail. 

As stated in Section 1, the mass and chemical composition of stars vary, 
and these are the characteristics which determine the history and lifetime of a 
star. New stars are continually born in a spiral galaxy, but they are not born 
from other stars, and consequently the mass of a new star is not in any sense 
inherited. The first generation of stars must have formed from clouds of almost 
pure hydrogen and helium. Subsequent generations of stars are observed to form 
from the cold, dense Giant Molecular Clouds of the interstellar medium. Given 
an environment which already contains surrounding stars, GMCs are only able 
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to remain cool because they contain dust, which acts as a shield to starlight, and 
because they contain organic molecules, which are able to radiate excess heat 
(Smolin pi 10). Organic molecules contain carbon, a by-product of the nuclear 
fusion processes in previous generations of stars. 

The star-formation processes within GMCs are thought to be triggered by 
the shock-waves from the supernovae explosions of nearby high-mass stars, 
formed within neighbouring GMCs 10 million years previously, or less. The 
shock-waves purportedly compress the interstellar medium, and instigate star 
formation. Because one spell of star-formation is triggered by the death of high- 
mass stars from a previous spell of star-formation, this process is referred to as 
self-propagating star-formation (Smolin pl28). Indeed, the inner lanes of spiral 
arms are often observed to be dark and dusty, while the outer sides contain the 
star-forming regions. This suggests the star formation is directed against the 
direction of rotation, which would explain why spiral arms rotate slightly slower 
than a host galaxy. 

Nicolson (p208) suggests that after a burst of such self-propagating star for- 
mation, differential rotation stretches the region into a spiral arm. He suggests 
that in a flocculent spiral, random bursts of star formation produce numerous 
spiral arms by this mechanism, each of which fades away after star formation has 
ceased. Subsequent star formation elsewhere in the galaxy then creates more 
spiral arms. This differs somewhat from Smolin's notion of self-propagating star 
formation, in which he states that "the waves of star formation neither die out, 
nor grow uncontrollably, but propagate at exactly the right rate to persist in 
the galaxy indefinitely," (1997 pl35). 

The debate between the two different theories which attempt to explain 
spiral structure is particularly interesting when one appreciates that they take 
opposite stances on the basic cause and effect relationships which operate in 
a spiral galaxy. The density wave theory represents a spiral to be a spiral of 
density waves, and explains star formation as the consequence of the density 
waves. In contrast, the theory of self-propagating star formation represents 
spiral structure to be the consequence of waves of star formation, (Tayler 1993, 
pl45). 

It would not be correct to say that the chemical composition of a new star is 

independent of the chemical composition of previous generations of stars. The 
chemical composition of the interstellar medium in a spiral galaxy is changing 
with the passage of time. The material expelled from one generation of stars 
provides a high-metallicity contribution to the medium from which the next 
generation of stars is composed. High metallicity purportedly inhibits the for- 
mation of higher-mass stars, hence the relative birthrate of lower-mass stars in 
a spiral galaxy increases with the passage of time. Accordingly, the birthrate 
of stars in the solar neighbourhood is claimed to be inversely proportional to 
mass m according to the Salpeter form, m~'^/^, (Tayler 1993, pl49). It is also 
suggested that the formation rate of high-mass stars is less than the formation 
rate of low-mass stars because the energy from a newly-created high-mass star 
heats up the medium from which it was born to such an extent that the star 
formation process halts, (Smolin pl27). In general statistical terms, the number 
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density Nm of type-m objects in such a population is the product iV„i = _B„i • Im 
of the type-TO birthrate Bm with the type-m hfe-time Im- If we let Nm denote 
the number density of mass-m stars, then it follows that because small mass 
stars have a greater lifetime, even if they have the same birth-rate as high-mass 
stars, they will come to dominate the population. The eventual domination of 
lower-mass stars is therefore the consequence of both statistics and the physics 
of star formation processes. In a spiral galaxy, however, where the birthrates 
are far from being spatially homogeneous, and the majority of star formation 
occurs in the spiral arms, one cannot infer birthrates from the known lifetimes 
of stars and the observed number densities in regions which are remote from the 
star formation in spiral arms. 

The population of stars in a spiral galaxy is therefore a type of population in 
which: (i) there are variable characteristics distributed within the population, 
and a small subset of these characteristics define different types within the 
population; (ii) each member of the population has a finite lifetime, determined 
by the values of the type-defining characteristics; and (iii) new members of the 
population are not reproduced from existing members of the population, i.e., 
new members of the population are born without having any parent(s) in the 
existing population. Suppose in addition that (iv) population members of each 
type are created at approximately the same rate, or short-lifetime members are 
created at a lower rate than long-lifetime members. Such a population evolves to 
be dominated by the long-lifetime objects. Such a population is neither evolving 
randomly, nor is it evolving by natural selection. 



References 

[1] Cox, J. P., Salpeter, E.E. (1964). Equilibrium Models for Helium-Burning 
Stars. III. Semi-Degenerate Stars of Small Mass, Astrophysical Journal, vol. 
140, p.485. 



[2] Disney, M.J. (2000). The case against cosmology, arXiv:astro-ph/0009020 
vl 1 Sep 2000. 

[3] Gabriel, M., Noels-Grotsch, A. (1968). Stabilit sculaire des toiles de carbone 
pur, Annales d'Astrophysique, vol. 31, p. 167. 

[4] Hansen, C.J. (1978). Secular stability: applications to stellar structure and 
evolution. Annual review of astronomy and astrophysics, volume 16, ppl5- 
32. 

[5] Hansen, C.J., Kawaler, S.D., Trimble, V. (2004). Stellar Interiors: Physical 
Principles, Structure and Evolution, 2nd Edition. Berlin-Heidelberg-New 
York: Springer Verlag. 

[6] Kahler, H. (1978). The Vogt- Russell theorem, and new results on an old 
problem, pp303-311 in The HR diagram - The 100th anniversary of Henry 
Norris Russell, A.G. Davis Philip and D.S.Hayes (eds.). Dordrecht: Reidel. 



16 



[7] Karttunnen, H., Kroger, P., Oja, H., Poutanen, M., Donner, K.J. (2003). 
Fundamental astronomy, Fourth Edition. Berlin-Heidelberg-New York: 
Springer Verlag. 

[8] Kippenhahn, R., Weigert, A. (1990). Stellar Structure and Evolution. 
Berlin-Heidelberg-New York: Springer Verlag. 

[9] Nicolson, I. (1999). Unfolding our universe. Cambridge: Cambridge Uni- 
versity Press. 

[10] Pasha, I.I. (2004). Density- wave spiral theories in the 1940s, I. 
|arXiv:astro-ph/0406142 

[11] Rindler, W. (2001). Relativity: Special, General, and Cosmological. Oxford: 
Oxford University Press. 

[12] Russell, H.N., Dugan, R.S., Stewart, J.Q. (1927). Astronomy Vol. 2. 
Boston: Ginn and Co. 

[13] Smolin, L. (1997a). The Life of the Cosmos. London: Weidenfeld and Nicol- 
son. 

[14] Tayler, R.J. (1993). Galaxies: structure and evolution. Cambridge: Cam- 
bridge University Press. 

[15] Tayler, R.J. (1994). The Stars: their structure and evolution. Cambridge: 
Cambridge University Press. 

[16] Vogt, H. (1926). Die Beziehung zwischen den Massen und den absoluten 
Leuchtkraften der Sterne, Astronomische Nachrichten, volume 226, p. 301. 



17 



